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Abstract 

A Ul(3) <8> Ur(3) low-energy effective lagrangian for the nonet of pseudogoldstone bosons that 
appear in the large N c limit of QCD is presented including terms up to four derivatives and explicit 
symmetry breaking terms up to quadratic in the quark masses. The one-loop renormalization 
of the couplings is worked out using the heat-kernel technique and dimensional renormalization. 
The calculation is carried through for Ui(ni) <8> Un{ni), thus allowing for a generic number ni of 
light quark flavours. The crucial advantages that the expansion in powers of 1/N C bring about are 
discussed. Special emphasis is put in pointing out what features are at variance with the SUl®SUr 
results when the singlet rf is included in the theory. 
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1 The U L (3) <g) U R (S) symmetry. 



The pattern of the lowest-lying states in the spectrum of strong interactions uncovers an approxi- 
mate continuous symmetry of nature, the so-called chiral symmetry, which is spontaneously broken. 
The octet of pseudoscalar particles - tt, K, and r\ -, with masses much smaller than those of the next 
excited states - the octet of vector particles p, to and K* , the baryons-, are the accepted candidates 
for pseudo-goldstone bosons associated to the spontaneous breaking of the symmetry. 

This approximate symmetry is well incorporated in QCD as three of the quarks happen to be 
light. In the (chiral) limit of vanishing m u , m^, m s the QCD lagrangian has the symmetry freedom 
of arbitrarily rotating with unitary matrices the quark field components in the space of flavours 
(u,d,s), independently in the Left and the Right sectors of chirality eigenstates. The symmetry 
group is Ul(3) (8) L^i?(3) and is explicitly broken by the light quark masses; if it had not, the lightest 
mesons would indeed have been massless particles. This breaking is small, though, since the light 
quark masses are much smaller than the typical hadronic scale of a few hundred MeV. (This is 
certainly so for the u and d quarks ( 2 < m u < 8 and 5 < < 15 MeV) and still approximately 
verified for the heavier s-quark ( 100 < m s < 300 MeV), Q). 

Empirically, however, only the SUl(3) ® SUr(3) symmetry subgroup, spontaneously broken to 
SUl+r(3), is manifest: instead of a nonet of light pseudoscalars only an octet is observed. Of the 
remaining Ul(1) <8> Ur(1), the vector part provides the conserved baryon number current, whereas 
the axial Ua(1) does not seem reflected at all in the spectrum, either as a conserved quantum 
number or as a goldstone boson. The first possibility would imply that all massive hadrons would 
appear in parity doublets and this is not what is observed. On the other hand, by its quantum 
numbers the rj would be the ninth candidate for goldstone boson; but if the Ua(1) were realized 
in the Goldstone mode and explicitly broken only by the same quark mass terms that break the 
SUl{3) ® SUr(3) one would expect that the ninth pseudo-goldstone boson would have a mass 
similar to the pion (to the masses in the octet): actually, a singlet pseudoscalar meson ought to 
exist with a mass smaller than y^m-Tr 0. Such a particle is missing in the spectrum. There is the 
rj instead but it is indeed so much heavier than the tt, K, rj that, at any rate, it seems hard to 
conceive it on the same footing as the light pseudoscalars, altogether in a nonet. 

This puzzle is part of what is known as the Ua(X) problem of QCD ||, and originates in that the 
Ua(X) symmetry of the QCD lagrangian is anomalous at the quantum level. Yet a conserved ninth 
singlet axial current can still be defined in spite of the ?7a(1) anomaly, it is not gauge invariant and 
't Hooft showed how its conservation could be broken by non-perturbative effects H|. 

Nevertheless, there is a limit of QCD in which the rj appears as the ninth, genuine goldstone 
boson: the limit of large number of colours N c . 't Hooft proposed a systematic expansion of 
QCD § @ with the inverse number of colours, 1/N C , as the expansion parameter. With the only 
assumption that in the framework of this 1/N C expansion QCD confines, a few qualitative features 
of the strong interactions already emerge by keeping only the leading terms. Of special interest to us 
is the result that the pattern of chiral symmetry breaking is exactly Ul(jii) ®Un{ni) — ► Ui+nini), 
where n; stands for a generic number of quark flavours @ , which resembles very much the pattern 
for chiral symmetry breaking observed in nature. Now in the N c — > oo limit there is no spoilt £7,4(1) 
anymore, since the anomaly in the divergence of the singlet axial current is 1/N C suppressed and 
the full Ul ® Ur is recovered along with the entire nonet of goldstone bosons |§ ||. 

The present study is based on the framework that the 1/N C expansion provides. Following the 
steps pioneered by Weinberg [ 10 1 and Gasser and Leutwyler [jnf |, we write the low-energy effective 
lagrangian of the N c — ► oo limit of QCD, with m u = = m s = 0: it is a chiral lagrangian that 
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involves the whole nonet of pseudo-goldstone bosons and is invariant under Ul(3) ® Ur(S). The 
departures from this scenario, which stem from the explicit breaking of chiral symmetry by quark 
masses and from the ?7a(1) anomaly, are treated perturbatively, in powers of the quark masses and 
1/N C . It is conceivable that a good picture of the lightest hadrons and their interactions at low 
energies could emerge from this approach. 

Many authors jl2j have discussed how to construct such lagrangian, i.e., how to extend the 
symmetry to ?7l(3) ® Ur(3) and properly take into account the effects of the £7^(1) anomaly at 
the same time, nicely organized in powers of 1/N C . In these articles the physical consequences to 
lowest orders in 1/N C and the derivative expansion have been worked out as well. We closely follow 
their work. 

The present study is devoted to extend the analysis to a full 0(p 4 ) chiral lagrangian, i.e., 
with terms kept up to four derivatives and quadratic in the quark masses. The conservative 
bookkeeping of quark masses as two chiral powers, m q = 0(p 2 ), is adopted, as in [11] (see also 
[14]). The celebrated Gell-Mann Okubo relations amongst the light pseudoscalar masses squared 
follow from this chiral power counting in a rather natural way. The exact mechanism of chiral 
symmetry breaking being unknown, however, many different possibilities other than m q = 0(p 2 ) 
have not been ruled out hitherto, either from the experiment or from QCD. None of them has 
been considered here. They can be adequately treated in the framework of Generalized Chiral 
Perturbation Theory as proposed in |15[ |, if needed. In order to discriminate amongst the various 
possibilities, the nature of chiral symmetry breaking needs to be established on more solid grounds. 

Our study is to all orders in the large- N c expansion - in a sense that will be later qualified. 

We calculate all one-loop divergences to the effective action using the heat-kernel technique 
and dimensional regularisation, and carry out to this approximation the renormalization of the 
couplings. 

The article is organized as follows: in section 2 the method of external sources and the generating 
functional are briefly reviewed and the notation is set. In section 3 the chiral lagrangian including 
terms up to 0(p 4 ) and the one-loop effective action are put forward, before the 1/N C expansion is 
performed. Next, in section 4, the calculation is organized in powers of 1/N C and the first non-trivial 
terms are also given. The conclusions and the appendices follow. 



2 The method of the external sources. 

In this section we briefly review the symmetries of QCD that are relevant for the chiral lagrangian 
in the large- N c limit. 

The QCD lagrangian can always be written with a diagonal quark mass term, as 

£qcd = ~\tt (C^GH + Y^fljf - m f ) Qf , (1) 

where / labels the light quarks qf that appear in ni number of light flavours. In nature, n\ = 3 
at most, for the flavours u, d, s. Although explicitly omitted, the quark fields also carry a colour 
index q c which labels the fundamental representation of the gauge group SU(N C ) of colour. Nature 
has chosen N c = 3. The covariant derivative acts on colour indices through the gluon matrix, 
Gfj, which form an adjoint (N 2 — l)-dimensional representation of SU(N C ). On the quark fields it 
acts in the usual way, diagonal in flavour indices, 

D,q c = d,q c -i^={G,) c cl q c ' . (2) 
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The field strength matrix is G^ = d^Gy — d u G^ — [Gfi, G v ], whereas G M = G^ ', the 
sum over a runs from 1 to N% — 1 and (A a ) stands for the N c x N c Gell-Mann matrices of SU(N C ). 

The heavy quarks c, b, t are also omitted from the lagrangian in ([!]) for they decouple from the 
strong low energy processes which only involve light pseudoscalars with u, d, s quantum numbers. 

The N c dependence that accompanies the coupling constant, is explicitly displayed. Apart 
from the usual combinatorial factors that appear in the Feynman diagrams, this extra dependence 
in iV c must be added in order to allow for a smooth, non-trivial iV c — ► oo limit of QCD. We shall 
revert to the issue of counting the leading powers of iV c of each effective coupling in section 4. 

If the light quark masses are switched off, the lagrangian (p]) becomes invariant under the 
symmetry group Ul{ui) <8> Ur(tii), n\ = 3. Henceforth we keep the number of light flavours, ni, 
generic in the expressions. This is most easily seen by writing the lagrangian above in terms of the 
quark left and right components, 

1 -75 1 + 75 f o\ 

QL = q , QR = q- (3) 

The terms in (|lj) that involve the quark fields read, 

q L (il^)q L + qR^D^qa . (4) 

The symmetry of rotating independently the components and ((/_r)/ in flavour space with 

unitary matrices is manifest. It is a global symmetry that is explicitly broken by the quark masses. 
However, if all the quark masses were the same there would still be an invariance under the diagonal 
vector subgroup UL+nini). In that case the subgroup coincides with the unbroken subgroup after 
spontaneous breaking of Ul(tii) <8> Un{ni). 



As it is well known not all the symmetries of the classical action are maintained at the quantum 
level; the quantum theory thus generates anomalous contributions to the divergences of some 
currents - they are no longer conserved. 

The low-energy effective action is a convenient bookkeeping device to encode the symmetries 
of the underlying theory - QCD -, which automatically incorporates all the unitarity features of 
quantum field theory. In writing the effective lagrangian care must be taken that all the (chiral) 
Ward Identities (WI) among the Green's functions are well implemented, including the anomalous 



ones. Actually, the method put forward in ref. |TT| constructs a solution to these WI's. It is based 



upon the transformation properties of the generating functional. 

The probability amplitude of transition from the vacuum in the remote past to the vacuum in 
the far future, in the presence of terms in the lagrangian that couple the external sources linearly to 
the currents, as in (g) below, contains all the Green's functions among these currents. Its logarithm 
Z[f] is the generating functional of the connected Green's functions, 



etZ[f] = I d Xl dx 2 ...dx n f£( Xl )f?^x 2 )...f£(x n )(0\T J&(si)j£( a *)..j£(s„)|0) 1 (5) 

11 Til 



J's and f's stand for generic currents and external sources, respectively. 

We shall consider both bilinear quark operators (currents) and the topological charge operator 
coupled to external sources, added to the QCD lagrangian, 



£ = £-QCD + qhl t J- tI '{x)q L + qR"l^(x)q R - qR{s{x) + ip(x))q L 

g 2 0± 

16vr 2 N, 



^(.s(,-)-ip(x)) 9fl -^^rr(G^G^), (6) 
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QV>v — e^ ua ^G a p. The first two new terms correspond to sources for the Ul ® Ur Noether currents 
of QCD, the non-singlets are the generators of Current Algebra; s is a source for the quark mass 
term and p for pseudoscalar bilinears with the quantum numbers of it, K, rj and rf. The sources 
I at r ui s an d p are hermitian ni x ni matrices; 9 is a real function. The axial and the vector Vn 
sources are defined so that L = — and r^ = v^ + a lJL . One can, formally, write the generating 
functional path integral, 

exp{iZ[l,r,s,p,9]} = J [dqdqdG^] exp{i J dxC}. (7) 

The connected Green's functions are obtained by performing functional derivatives of Z with respect 
to the sources. 

In order to further constrain the form of the effective lagrangian it is a convenient trick to 
promote the global Ul 8) Ur transformations - that leave the QCD lagrangian invariant - to local 
ones by allowing the external sources to transform along with the dynamical fields. The combined 
set of local Ul <8> Ur transformations^], 

g L (x) G U L , 9r{x) G Ur 

QL -> 9LQL 

QR -> 9RQR 

Ifj, -> 9l1^9 l + igLd^gl 
Tfi + WRd^gR 

s + ip -> gij (s + ip) , (8) 

naively becomes a local gauge symmetry for the lagrangian density in eq. (1). This is not quite so, 
due to the fact that the transformations in (||) also induce a non-trivial anomalous Ua(1) transfor- 
mation on the generating functional Z[l, r, s, p, 9]. Its origin may be traced to the transformation 
properties of the fermionic integration measure |l7j (or, if one wishes, of the fermion determinant) 
once it is properly regularized. This is reflected in the anomalous divergence in the ninth axial 
singlet current, 



16vr 2 N c 

This anomaly-related effect turns off any potential advantage inherent to the existence of a gauge 
symmetry which eventually would severely constraint the form of the effective action. However, 
this drawback is obviated as the J7a(1) anomaly contribution may be altogether eliminated by 
judiciously choosing the transformation law for the external field 9(x). Indeed, for infinitesimal 
gL = I + i(f3 — a), gR = I + i(f3 + a), the source 9{x) ought to change as 

9{x) ->0(x)-2{a(x)), 

(here we have switched to the standard notation and denote the trace operation over flavour indices 
by brackets trp(...) = (■■■)) in order for the £^(1) anomaly to cancel. The term generated by the 
anomaly in the fermion determinant is explicitly compensated by the shift in the 6 source. 

A subtlety is still to be analized. The set of local gauge transformations we have constructed 
in (||) also induces a non-abelian anomaly. This new drawback can not be circumvented and 

^Being unitary, gn and gL can be always parametrized as ga = exp(if3)exp(ia) , gL = exp(if3)exp(—ia) , with a 
and /3 ni x ni hermitian matrices. A pure vector transformation has a = 0, whereas an axial one has /3 = 0. 
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needs explicit consideration. As discussed in [16], imposing upon regularization the requirement of 
conservation for the nine vector currents, the change in Z under (|8|) reads 



5Z = — J dx (a(x)Q(x)}, 

4 2 8 4" 

VapVpv + g (V a a |fl )(V^a I/ ) + -i{v a p,a,j,a„} + -la^v^ay + +-a a a / 3a /1 ai, > 

w a/ 3 = o>a«/3 - dpVa -i[v a , vp] , = d a ap - dpa a -i[v a , ap]. (10) 

The terms in 5Z appear due to triangle AVV Feynman diagrams which involve insertions of an 
axial and two vector quark bilinear operators. Higher polygon-shaped diagrams, quadrangles and 
pentagons are anomalous as well and give rise to the cubic and quartic terms in (|l0|). Adler and 
Bardeen showed that the coefficients of the anomaly are not affected by higher-order radiative 
corrections, i.e., diagrams with more than one-loop do not contribute to the anomalous terms |jTq 1. 
5Z fulfills the Wess-Zumino consistency conditions fl9[| . 

The integrated form of this anomaly must be added by hand to the low-energy effective field 
theory that we shall construct. Such a theory will consist of interacting bosons and therefore will 
not be anomalous. Thus, the above effect will be contained in an additional term to the effective 
lagrangian. 

It is worth pointing out that unlike for SUl <S> SUr, here there is the possibility of combining 
one Wess-Zumino- Witten (WZW) vertex (which start out at 0(p 4 )) with the 0(p°) pieces in the 
lagrangian and generate additional one-loop divergences at 0(p 4 ). They will be given elsewhere 
[p0| . Notice that for SUl®SUr this cannot occur at 0{p A ) because the lowest chiral power counting 
are 0(p 2 ) terms, which generate divergences at 0(p e ). 
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3 The chiral lagrangian 

The mere knowledge of the symmetries and the way they are realized provide an enormous insight 
for they are reflected in every aspect of the theory, e.g., in the spectrum, the Green's functions, 
the interactions, to mention a few. In the case of the strong interactions the symmetries are the 
bulk of the information which is available at low energies, since QCD is non-perturbative there. 
Although a lot of effort has been put in numerical calculation projects and enormous progress in 
solving the technical difficulties has been achieved, so far the problem has remained too hard to 
tackle satisfactorily. 

In this section we shall write an effective lagrangian for the soft interactions of the lightest 
particles in the spectrum, the nonet of pseudoscalar mesons, tt, K, 77, rj . 

The effective lagrangian is a combined statement about the degrees of freedom and the symme- 
tries that are relevant for the processes under study. Effective refers to the choice of field variables, 
in this case fields for the tt, K, 77, rf particles that are observed in the range of energies below 
the p-meson mass m p . More generally, being the lightest particles in the spectrum the use of the 
effective lagrangian will be the determination of the long distance behaviour of any of the QCD 
Green's functions, where they are expected to dominate. 

Being a symmetry statement, the strategy consists of writing down for the effective lagrangian 
the most general expression that contains all the independent terms compatible with the symme- 
tries, multiplied by unknown constants. The fate of many effective theories is to be of little practical 
use if the number of unknown constants blows up. They can always be fitted from experiment but 
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too often the number of experimental results available is of the same order (if not smaller) as that 
of constants, rendering the approach little predictive. 

In the present case the spontaneously broken symmetry character and the consequent goldstone 
nature of the it, K, rj, rj impose severe constraints on the form of the interactions. The number 
of unknown constants reduces in a drastic manner if one restricts to the lower terms. Actually, 
they would reduce to a handful of them if there were no C/a(1) anomaly effects to incorporate, as 
happens in the SUl <8> SUr lagrangian if one only keeps terms up to 0(p 4 ). Fortunately, the swarm 
of new terms that the Ua(X) anomaly introduces carry high powers of 1/N C , and to lower orders in 
\/N c only a few survive. 

Following [11 1, [Fyl we collect the nine pseudoscalar fields in a hermitian matrix 3>(x), 



<J>(x) = r]°(x)\ + ir(x) 

where tt(x) = tt(x) • A, Ao = \J~^I an d A are the Gell-Mann matrices of SU(ni). For ni = 3 (see 
Appendix B), 



<J?(x) 



( ^ + + Jlr? n+ _ K+ 



(11) 



The unitary n\ x n\ matrix U(x) is the exponential of &(x), 

U(x) = e^ x)/s , (12) 

/ is an order parameter of chiral symmetry breaking and gives the strength of the coupling between 
the goldstone bosons and the currents that are spontaneously broken and do not annihilate the 
vacuum. In this case det U is a phase, 

det U = expiiV'Mi rj°/f)- 
Under Ul <S) Ur, U transforms linearly, 

U g R Ug[. (13) 



The transformations induced by ( [l^ ) on $ are more involved. Under Ur + l(1), $ is automati- 
cally invariant, as it should, for mesons must carry baryon number equal to zero. Under SUl+r(3), 
transforms linearly: it contains two irreducible representations, the octet ir and the singlet 77°. 
Under an axial transformation 

U — » e ia Ue ia , 

$ changes nonlinearly, 

$ + 2a + 0(a 2 ). 

This can be understood on geometrical grounds since the fields in $ may be regarded as coordinates 
that span the coset space Ul®Ur / Ur + l, upon which the Ul®Ur acts: the fields themselves are, 
in a sense, parameters of a group element, and the nonlinearity reflects the group transformation 
law when written in terms of the continuous parameters. For (log U) = i\/2nirj° / f , 



(log CO — » <logl7> + (logC0itfi)> = (logtf) + 2i< 
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rj gets thus shifted only under an axial transformation and is invariant under any vector one. 
Under Ul <8> Ur it never mixes with any of the it components. 

With this choice of fields, the origin of the derivative couplings among goldstone-boson becomes 
transparent: a lagrangian invariant under global Ul (8> Ur ought to reduce to zero if U were a 
constant matrix, for then, by virtue of the symmetry, it could always be transformed away with a 
global unitary rotation: the couplings need thus be derivative couplings j^J. The expansion will 
be in powers of momenta of the soft light mesons divided by a scale of chiral symmetry breaking, 
which is of order ~ 4irf ~ m p , the mass of the next excited state, the p- meson. In addition, in the 
present case the Ua(1) anomaly introduces novel couplings of the tjq meson that are not derivative 
couplings. 

The objects at hand to construct the chiral effective lagrangian are the matrix U and the 
external sources, r^, In, (s + ip) and 9. It is useful to introduce covariant derivatives for the fields 
and the sources. The covariant derivative will act on flavour space and, as usual, its action will 
depend on the transformation law of the object it derives, 



D^U 



d»U 



D^{s + ip) 
D^logU) 



ir^U + iUlfj,, 
dfj,(s + ip) - ir^(s + ip) + i(s + ip)^, 
(U^D f ,U) = d tl (\ogU)-i(r tl -l^, 



One is led to introduce field strengths for the vector and axial sources 



^R 



Under local Ul <8> Ur they transform as 

D^U 
D^(s + ip) 
F L 

F R 



du,r u - d v rn - i[rn,r v ] 



(14) 



(15) 



9R (D^U)g[ 



(D^s + ip))g L 



9R 



9lF^ 9l 
9R F ?u9r- 



The combination 



X(x) = (logU(x)) +9(x) 



■ V^nj o , 

' — " + 



= it 



(16) 



(17) 



is invariant, and so is any function of X [13|. This is a novelty of Ul <8> Ur and it is possible 
because (log U) does not vanish, as it does for SUl®SUr. Due to the U^iX) anomaly, each invariant 
operator generates, in reality, an infinite set of invariant operators, since the symmetry allows to 
multiply it by any function of X and still remain invariant. Notice that this method of finding new 
operators by multiplying the old ones by functions of X never introduces new derivatives to the 
vertices. It is for the same reason that counting the number of derivatives and retaining operators 
that contain up to a certain number does not limit the number free constants, as it used to in 
SUl ® SUr. At each order in the derivative expansion we find an infinite number of constants. 
It is customary to introduce the source x( x ) 



X = 2B(s + ip), 
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which transforms as the U matrix. B is a constant that is related to the quark condensate (qq). 
Its relevance comes from the fact that the explicit symmetry breaking driven by the quark masses 
provides the former goldstone bosons also with a mass. If m q denotes a quark mass, to lowest chiral 
order there is a contribution to the meson masses which involves B and is given by m 2 = 2Bm q . 
We assume that it is the bulk of it: we make the hypothesis that further pieces which involve order 
parameters other than B are negligible, comparatively. This leads to counting x as 0(p 2 ). The rf 
gets an additional piece to its mass through the ?7a(1) anomaly, which is 0(1/N C ), and which does 
not add to the masses in the octet (^). The singlet-octet mass splitting is very big in practice, of 
~ 400 MeV in the least extreme case. 

All the constants that appear in the effective theory are free, to be fitted from experiment. 
Although we have not been able to compute them from QCD, there are exact inequalities that 
have to be verified, which are based on the vector structure of QCD. These relations are non- 
perturbative |^q] . 

From the effective lagrangian, one can make contact with QCD through the generating func- 
tional Z[l,r,s,p,9] introduced in section 2. One demands that the same functional - the same 
Green's functions - should be obtained by starting from the effective theory as from QCD. It can 
be formally written in terms of the light pseudoscalar fields, collected in U(x), as 



exp{iZ[l,r lS ,p,e}} = I [dU] e i S** L '" 



low modes 



(18) 



The chiral lagrangian only copes with the long distance behaviour of Green's functions. Only the 
lowest modes have physical significance in (|l~8|). For distances smaller than l/m p the approach is 
inappropriate and the integrals over loop momenta have a natural cutoff associated with them that 
The higher modes, corresponding to more energetic pseudoscalar s, and the rest of massive 



is m 



hadrons are integrated out and their effect manifests through the coupling constants of the effective 
theory §7|. 

The most general lagrangian invariant under (local) Ul (g> Ur [] that includes terms with two 
derivatives or less, or one power of %i is the following [pjj , fl~3|| , 



C 



(0+2) 



- W (X) + W\{X){D I JJ^ D^U) + W 2 {X){U ] x + X ] U) + iW 3 (X)(tf X ~ X ] U) 
+ W 4 (X)(tfD^U)(tfD»U} + W 5 (X)(tf(D^U)}D»§ + W 6 (X)D tl §Di*§. (19) 



Parity conservation implies that they are all even functions of X except for W3 that is odd. 
Furthermore, VF4(0) = 0, Wi(0) = VF2(0) = 4 gives the correct normalization for the quadratic 
terms. 

The first term Wq has neither derivatives nor powers of x, and therefore counts as 0(p ). The 
rest of the terms count as 0(p 2 ). The 1/N C power counting of the diverse couplings is given in the 
next section. The term Wq brings a mass terms for the rp 



m 2 v0 



-^w?(o), 



(20) 



2nii 

UaW ' f 2 

"Terms of the sort '^f a _ Q {\ a U^ D tl U)(\ a U^ D^U) and alike, given the properties of the {A a } matrices, can be re- 
expressed in terms of the operators written in the text. In this particular case, with the relation y^_ A^A£j = 25uSjk, 
it becomes -2((D M [/ t )(D"[/)) 
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whereas the first term in the expansion of W$ gives a contribution to singlet-octet mixing from 

Ua(1). ' 

The N c — > oo limit of QCD actually imposes more restrictions on the effective lagrangian than 
the symmetry Ul®Ur alone. Under Ul®Ur, the fields rp and tt never mix their components. There 
is no reason why the particles created by them should bear any sort of relationship whatsoever, 
as though they belonged to the same irreducible representation, like the ir do. There is no reason, 
either, why in the definition of U(x) in ( |i~l| ) rp and 7r should appear in the exponent divided by 
the same constant /, i.e., normalised in the same way. One could have written instead 

*(\ 

U(x) = e\V ntf0 s > , 

with / and fo unrelated and, yet, ( |i~9| ) would be Ui (g> Ur invariant. The proper way to cast this 
issue requires to fix the normalization of each field by looking at the kinetic energy too, not at U{x) 
only. In (|19|), the kinetic energy terms are ^(D^D^U) + W A (0)(U ] DJJ)(U ] D^U), which read, 

2 (W • ( a n ) + 2 ( /2 ^ J d ^od^o- 

The normalization condition (/ 2 — AniW^O)) / f$ = 1 relates three constants; it may be viewed 
as an arbitrariness in their definition, for a change in /o can be always compensated for by an 
appropriate change in W4(0). Once this normalization is fixed, the strength that the singlet rp 
couples to the singlet axial current is /o whereas the octet couple to the octet axial current with 
strength /. 

Unlike Ul ® Ur, which does not have a dimension nine irreducible representation, the N c — > 00 
really enforces a nonet symmetry, with the ir, K, r/, rf all having identical properties. Indeed, the 
planar Feynman diagrams that contribute to (jf\{x)jfl(Q))) and to ( J^(x) J^J(O)) (a = 1, 8), 
with Jg = qjn^^-q, are the same since the qq — > gluon annihilation diagrams that would only 
contribute to the singlet channel turn out to be 1/N C suppressed (OZI violating processes). Barring 
for 1/N C corrections the two decay constants coincide f / fo = 1. Moreover, if the same quark mass 
were switched on for all light quark species, a mass would be generated identical for the singlet 770 
as for the octet particles || Q. 

The standard normalization of the kinetic energy and nonet symmetry require W4(0) = 0. 
Without loss of generality, in that case the term W^(X) may be eliminated altogether by a change 
of field variables of the type U — > U exp[iF(X)], that maintains the transformation properties for 
U (but changes the normalization conditions of rp). 



The methods of the background field and the steepest descent, when applied to the functional 
integral (|l8l), provide the loop-wise expansion of the generating functional. The Green's functions 
are read off from the generating functional, about its minimum when the external sources are 
switched off, x = 25 diag(m u , raj, m s ) and 9 = 0. To lowest order we assume that the minimum 
is achieved at Uq = I, which is compatible with the equation that minimizes the lowest order 



effective action jll|. 



In order to include one- loop corrections, one proceeds to introduce a background field matrix 
U(x), and expand ( p~9| ) about this background configuration. For that we decompose U(x) as 

U(x) = U(x)S(x), E(x) = e iA(x) , 
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with A the matrix of quantum fields PI By choosing U to transform under Ul <8> Ur as U, the 
transformation laws for E and A become 

A -> ffLAg^, S -> ffiSg[, 
and, therefore, their covariant derivatives are 

= 0„£ - i[Z M , S], Z^A = 0„A - i[Z„, A]. 

Next we expand the lagrangian in (|l^) about U and keep terms up to quadratic in A. The 
one-loop effective action is obtained upon evaluation of the functional determinant of the dif- 
ferential operator that appears in the piece quadratic in A. Finally the background field U(x), 
which is a priori independent of the sources, is judiciously chosen so as to verify the equations 
6Z[U, I, r, s,p, 9] = with the sources held fixed. The method ensures that with this U, Z is the 
effective action T[U] - the generator of the one-particle irreducible Green's functions of fields gath- 
ered in U -, in the presence of the external sources. For the one-loop effective action it suffices to 
use the tree level equations of motion. The corrections to it would modify the two-loop effective 
action. 

Expanding the effective lagrangian in (|i~9l ) about U, disregarding the terms linear in A(x) that 
vanish with the equations of motion, it yields: 

£ + 2 (£/£) = C 0+2 (U) + A (A) 2 

+ 2W[(X)(A)(C<* D,A) + Wi(X)(D„A D^A) + W\ (X) (C^ [A, D^A]) 

+ W^(X)(A)(AN) - hv 2 (X)(A 2 M) 
+ iW^X)(A){AM) - ~W 3 (X){A 2 N) 

- W^X){A){D^A)D»0 + O{A)\ (21) 



where 



A ee \w^X) + \w^X){C^)- l -W^{X){M) 

^^'(X)(AT)-i<(X)(C M )W-i<(X)^W, (22) 



and 



ee Wd^u, m = W x + x j u, n = W x - x j U; 

all transform as C M , M -> ghMg\, N -> giNg^; also, C+ = -C^, = M and iV 1 " = -N. 
The functions Wi(X) have been Taylor expanded about the background value X = (log U) + 9 as 
follows, 

W fc ((log(tfE)) + 9) = W k (X) + iW> k {X){A) - I<(X)(A) 2 + 0(A) 3 . 



"The procedure is not manifestly left - right symmetric. This is not a worrisome issue given that the quantization 
of scalar fields in four dimensions does not have any anomaly that would favor one over the other. Simplicity reasons 
have lead to our choice. A more symmetric treatment would lead to the same final results. 
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Integrations by parts have been performed where necessary and the total divergences have been 
discarded. The equations of motion can be read off from the terms that are linear in the variation 
A, 

1 Wn 1 Wi Wi 

M 2 Wi 2 W% M Wi p ' 1 ^ ' 

1 Wo iW% 1 Wo i Wi 

+ Z_1tv + - ——M — ——^-(M) — - —^-(N) 

2 Wi 2 Wi 2 Wi V 7 2 Wi W 



+ 



1 Wr 1 Wfi \ - - IK 

2Wi 2 Wj V M A ; 2 Wi ^ V ; 



(Henceforth the arguments X are omitted from the functions W^'s and their derivatives that 
appear in the calculation; also, bars are suppressed from X and U). 

In order to be able to use the expressions collected in the Appendix A for the evaluation of the 
one-loop divengent parts it proves useful to perform a change of integration variables so as to leave 
the operator that acts on the quadratic piece in (|2l|) in the usual form, with the laplacian piece 



dud* multiplied by a constant, not by a function. For this purpose we change variables to 

A(x) = £^=<P(x), (24) 
and expand the hermitian matrix <p(x) in the basis of matrices X a (see Appendix B): 

Retaining the piece quadratic in the quantum fluctuating fields ip a , one finds 

Quadratic = ^ + ^ ^ (25) 



where 
and 



[d^ = ^+<V, (26) 



f = \{% + \c,)[\ a , \% + ^ ((C,\ a )(\?) - (C^)(X a )) . (27) 



Notice that <jf is antisymmetric in a, (3. It is this property what allows to integrate d a by parts 
as a whole, as though it were the single d u . For that reason the operator d u d* + a is manifestly 
hermitian. The expression in (Eq) differs from that of (|2l|), after the change of variables, by a total 



derivative. The evaluation of the Gaussian integral involves the expression (25), though: it is the 



determinant of the differential operator that is hermitian the one that has to be evaluated [n], 

t1 dn practice this means that any term of the sort ip a f^ 13 (x)d t± ifp (which, in general, is not hermitian as can be seen if 
one tries to bring the operator act on ip a , on the left) can always be written as <p a aff (x)d tJ '(p 13 — ^i/3 a 9 M (s^ /3 (a;)) ip 13 + 
^dfj. (y> a sff (x)^^ . Notice that the last term is a total derivative which we shall discard. The remaining part 
a r ^{x)d y ' is hermitian now. The aff{x) and sff{x) are the antisymmetric and symmetric parts of f^ l3 {x), respectively. 
This decomposition is unique. 
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The curvature associated to this connection is 

~ |(^) 2 (^A a )(aA^)-(^)(aA a )), (28) 
fj,, v are space-time indices, whereas a, (3, 7 label the Gell-Mann matrices of U(ni), 

QV» = + rfFgU - l -[C», C u ], 

and 



H " v = \m ~ Km)) imcu ~ {cncn + m " (m) ^ D " § " c ^ 



For a we find 



+ i((TA a )(A' 3 ) + (rA' 3 )(A Q )), (29) 



where 



1 fWi' 1 {W{\ 2 \ . , IE' 

s = -2 - 2 (wt) J + » 2 - iwr'^ + D » OT >• 

t = ^ + ^-2y) wc '-2 ( ' ; ' c " ), ^ ,B ' ,; ' ) 

1 Wo' Wl' Wo Wl Wl' 1 Wl 

1 (Wa Wl'\ - o 

R = w l M + l w l N - (30) 

The one-loop effective action is obtained by including the quadratic fluctuations about the 
configuration [U] that, consistently, minimizes the effective action itself. One needs to evaluate the 
integral of a Gaussian functional, with the known formal result 



J ^/det (d^ + a) 





which, upon exponentiation, contributes to the effective action as 

T One-loo P[tj] = J d A x £(0+2)( tf) + % _ Tr log + a) + J d A X C (i) {U). (31) 
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A word is needed on the proper definition of the previous expressions. The heat-kernel technique 
has been used to define the determinant (see Appendix A), and the divergences have been dealt 
with dimensional regularization. In order to absorb the infinities that result from the functional 
determinant, counter-terms of 0(p 2 ) and 0(p ), as well as new terms of 0(p 4 ) need be included. 

The determinant from the change of functional integration variables in ( |24T ) gives a contribution 
which is proportional to a singularity 5^(0), and dimensional regularization sets it equal to zero. 
(A similar remark was in order when a change of field variables allowed to cross out W^X) from 
the lagrangian in fllOD). 



We find, 



2 I 8 32 



{wd) (WW *) + ^(C^)(C U C») - \(C»C»C V ){C") 



+ I (^) 2 (WO - (C^C v cn) + ^(R) 2 + ^(R 2 ) + f S 2 + \{T) 2 
+ ^(r 2 ) + S(T) + | ((^) 2 - l) (WR) - \(C,C»){R) + \{C»){C»R) 



/<(""') + ~T'^{R) + rABT) + ( ^ I' W(C2% 



and 
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1 ^<„ MH ^ + driL\' 



+ T2{WJ {C >> Hn{C " ) ' 

which are the only structures that get divergent contributions at one- loop (see Appendix A). 

Since quantum scalar fields in four dimensions do not generate any anomaly to the Ul® Ur 
symmetry, the new terms needed to renormalize the one-loop result are necessarily in the list of all 
possible operators of 0(p 4 ) invariant under Ul <8> Ur. 

The list of independent operators is given below. The criteria used to select this particular set 
are the following: terms involving the derivatives D^M, D^N, D tJi F^ u and alike; three derivatives 
of § (D^D^Dj), (D^D U D I1 §); D^D^C Vl (D^Dj) or D^C can be removed as combinations of 
those in (^) plus terms with the piece D^C P : finally, these can be eliminated with the equations 
of motion (|23|). The two derivatives of 9 can always be chosen to appear under the form D^D^O. 
The rest of operators that are not independent have been removed upon integration by parts and 
with the help of the identities 

D n c v _ D » C v = _[ C r/* j c v \ + iff - irfFgU, 

[D»,D v ]C p = -i[F p \C p ], 
[D P ,D V ]6 = i(F R iU - F pv ). 
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All the following terms are real Oj = { . The first ones correspond to the twelve 0(p ) operators 
of SU L ® SUr (recall that = WD^U, M = W X + X f U, N = W X ~ X ] U and § = i0), 

O = (D^U D U U^ D»U D v tf) = {C»C V C^C V ), 

Ox = (Dprf D»U) 2 = {C»C^){C U C V ), 

2 = (Dptf D v U){D»tf D v U) = {C»C v ){C il C v ), 

3 = (D^U^ D^U D U U^ D U U) = (C^CpCC,,), 

4 = {D^U* DHU){rtx + X*U) = -(CC lt ){M), 

5 = (D^tf D»U(tf X + X f U)) = - (CC^M), 
Oe = (U^ X + X t U} 2 = (M) 2 , 

7 = (Uh-X^U) 2 = (N) 2 , 

8 = {x ] U X ] U + U^ X U ] x) = \{M 2 + N 2 ), 

O g = -i(F^ D^U D V U^ + Fl v DJJ* D V U) = i(C^C u (f^ + U^Fj^u)), 

O10 = (U^F^UFl^), 

On = (Fr^F^ + F L/JH/ F^), 

O12 = (xh) = \(M 2 -N 2 ), 

The following eight operators are obtained from the previous twelve by splitting up single traces 
into products of traces: (C M ) does not vanish, (C M ) / 0, for Ul <8> Ur. They read, 



Ol3 


= (C/ f L> M C/)(C/t D^U D U U ] D V U) = 


= -{cn(c^c u c u ), 


O14 


= (C/t DJJ)(tf Df"U)(D u U^ D U U) 


= -{C»){C lx ){C»C v ), 


Ol5 


= (C/t £> M C/)(C/t D V U){D^U ] D U U) 


= -(C)(C)<C M C,) ) 


Ol6 


= (C/t D*U){rf DJJ){tf D»U){tf D U U) = {C»){C^){C V ){C V ), 


On 


= (C/t D"U) (C/t DJJ) (C/t x + x t[/) 


= (^)(C M )(M), 


Ois 


= (C/t DJJ){D»tf X-D^U X ] ) = 


-(C^)(C M M), 


O19 


= (F R ^)(F^) + (F L ^)(Fn, 




o 20 


= (F R ^)(Fn. 





So far, all are parity even operators. The next seven are similar but have odd parity, 

021 = -i{D^ D»U {Uh-X ] U)) =i{NC»C^), 

022 = ~i(D^ D»U))(tf X -X f U) = i{N)(C»C lx ), 

023 = i(U* D^U){D^ x-D^U x t )=i(NC^)(C lx ), 

24 = D»U)(rf DpUKtfx ~ X j U) = z(iV)(C^)(C M ), 

25 = i{U^ X U ] x-X ] U X ] U)=i{NM), 

26 = i((Uh?-(x ] U) 2 ) =i(N)(M), 

027 = (C/t D,U)(Fl»rf D U U - Fjf A,C7 C/t) = (C^XC^Ff -tfFg'u)), 
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Three operators involve the t pupa tensor. The first two of them are odd under parity and O30 
is even. 

28 = e^F^tfFgU), 

29 = ie^ pa ((F£ u + tfF^u)c p C°), 

30 = W^l" -tfF£»u)c»)(C°), 

(32) 

From O31 to O57 they involve derivatives of the source 9, and are given in Appendix C. 
The operators that appear at 0(p 2 ) in (f2~l|) are 

E 1 = -(CC lt ) E 4 = (C' t )(C ft ) 

E 2 = (M) E 5 = (C»)D p 6 (33) 

£3 = t(JV) E 6 = D^eDj 

The effective lagrangian including up to one-loop corrections is, thus, 

6 57 

c One-ioo P = - W £(X, fj.) + W[(X, n)Ei + n)0, + finite non - local. (34) 

i=l i=0 

The following structure of counter-terms renders it finite (see Appendix A): 

5Wi(x) = nw l {l) (x,fi) + nn i (x^)x 00 + o{h 2 ), i = o,...,6, 

A(X) = (3l(X,ti) + hB i (X,ti)\ oo +0(h 2 ), i = 0,...,57. (35) 

with 

^^(i^ 1 *-^ 1 ')' <36> 

so that 

W[(X,fi) = W i (X) + hw[ 1) (X, f i), i = 0,...,6. (37) 

At this point we have reinserted the powers of Ti to help the counting of loops. Recall that the 
one- loop effective action carries one power of fi. 

The roster of functions JVs and B^s is given in Appendix D; they are the main result of our 
paper. The counter-terms in ( |35| ) are written in terms of the functions Cli's and B^s for the sake of 
concision. This notation, however, may seem a bit contrived. It should be read in the usual way of 
perturbation theory, namely with the functions JVs and B^s understood as their series in powers 
of X. The renormalization of the functions means the renormalization of the coupling constants 
which are the coefficients in these expansions. 

Parity, charge conjugation and time reversal ought to be conserved. Only the operators that are 
invariant under charge conjugation themselves can appear in the lagrangian. This is because X is 
invariant under C as well (see Appendix E). The list of C-violating operators is given in Appendix 
C. 

The result is valid for any value of n;. However, depending on the specific n\ considered, there 
are n^-dependent factorization relations among the traces of products of n\ x n\ matrices that are 
of relevance to us since some operators in the list become redundant, i.e., some can be written in 
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terms of a smaller subset. These relations follow from the Cayley-Hamilton theorem and have been 



extensively used in 11]. For ni = 3 they boil down to [21] 

A 3 - (A) A 2 + l - {^A) 2 - {A 2 )) A - det(A) = 0, 
for any 3x3 matrix, be it hermitian or not, and to 

£ (A X A 2 A Z A A )- <^Ma>(44>- E (^2X^3 At) 

6 perm, 8 ■perm 3 perm 

+ £ <^i4j> (4»> (A) - (Ai)(A 2 )(A3><^4) = 0. (38) 

6 perm 

The first relation ensures that the determinant of a matrix can always be expressed in terms of 
traces and justifies why the determinants of products of U matrices and their derivatives need not 
be considered independently. From the second relation, with 

A x = A 2 = C„, A 3 = A A = C v , (39) 

and summing over the indices [i and one gets 

2O -O x - 20 2 + 40 3 + 80 l3 - 20x4 - 40is - 16 = 0. (40) 

Thus, for m = 3 we can spare Oq and the 0{p i ) lagrangian will be written as 

57 

£ 4 = E L *( X ) *> ( 41 ) 

i=X 



where 



Lx 


= 01 


+ \Po 


Lx3 


= faa ~ 4/3o 




= 02 


+ 0o 


L14 


= 014 + 00 


L 3 


= 03 


-20o 


Lx5 


= 015 + 2/3 








Lx6 


= 016 + \Po 



(42) 



and Li = 0i for the rest. 



From Ox to Og the same notation as in SUl x SUr [11] has been kept, also for the coefficient 



functions. The constants Hx, #2 in [11] have turned into the functions Lxx{X), Lxi{X). The new 
operators that appear in this Ul (8> Ur lagrangian are labeled from O13 onwards and the coefficient 
functions Lj(X)'s follow suit. 



If one disregards all the coefficients associated to ?7a(1) one finds for the B^s in ( |35|) 

B o = l B i = I fi 2 = ^ B 3 = §| B A = \ B 5 = f B 6 = ± 

B 7 = 6 8 = a B 9 = a 6 10 = J Bxx = -^ 5 12 = | Bx3 = ± (43) 

Bx4 = Bx5 = Si6 = Bx7 = Bxs = — j Bxg = 24 -B20 = J2 

O30 involves an c^i/pQ- and does not need to be renormalized. 

In the case of ni = 3, taking into account that the same relations from (42) should hold, and 
renaming the constants as IYs, we obtain 

r i = i r 2 = re r 3 = ^ r 4 = g r 5 = § r 6 = ^ r 7 = o 

T8 = T6 ^9 = 4 Tio = —4 Tn = — g Ti2 = g Ti3 = Ti4 = jg (44) 

Ti5 = I r 16 = 32 r 17 = r 18 = — \ r 19 = ^ r 2 o = j2- 
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They coincide to those of SU L (3)®SU R (3) [|TTJ except for the terms that involve (M 2 ), (M) 2 , (N 2 ), 
(N) 2 : Tq, Ts, Ti2- The reason is that among the building blocks that have been used to write the 
chiral lagrangian, C M and F^ U R have vanishing traces in the case of SUl{3) (g) SUr(3), whereas 
neither M nor N do. Although it is less immediate to retrieve the SUl(3) (g> SUr(3) coefficients 



,[SU] 



11 

144 



JSC7] 



A. r 

48 1 12 



[SU] 



_5_ 

24' 



(45) 



from our result, there are simple relations that have to be verified. For instance, if one adds all the 
divergent pieces that go with the operators Oq, Os, Ou, sets x = m2 I f° r simplicity, and expands 
the operators, it is easy to check that the SUl(3) (g) SUr(3) and the Ul(3) <g) Ur(3) coefficients of 
(x X) an d v? 2 verify, respectively, 



( 



i 2 r[f u] + 2t[ su] + rtf i) 



:i2r 6 + 2r 8 + r 



12 J 



and similarly 
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[SU] 
6 



+ r [ 8 su] ) 



(3r 6 + r 8 



One recognizes the ratio | as the fraction of degrees of freedom in the two theories, for these 
coefficients multiply one-loop divergent pieces which in the two cases stem from tadpole diagrams, 
which give a constant divergent contribution for all the virtual mesons that travel around the loop. 
Therefore, the total result is proportional to the number of degrees of freedom that in each case can 
circulate. (Of course the same argument goes through for any number of flavours ni and a similar 



result holds in general. The SUiini) <8> SUr(ui) coefficients are 



,[517] 



2+nf 
16nf ■ 



r 



[SU] 



16n ; 



The second relation holds in the form n ' r6 -j~ rs 



„ r [SU] r [su\ 
niL R +1„ 



- . The first relation, that now involves 

the combination AuiTq + 21^ + Ti2, reduces to the previous one if one realizes that Ti2 = 2T% in 
either case.) 



One important difference between the SUl <8) SUr case and ours is that in the first theory the 
meson masses do not get any infinite contribution and in this case they do. This statement needs 
some qualification for the language it uses is the customary of renormalizable field theories, where 
the divergences that are generated require a fixed number of counter-terms, of same type as the 
terms in the lagrangian only. The chiral expansion in increasing number of derivatives is not a 
theory of this kind, rather it is non-renormalizable, because at each higher loop new terms are 
required to absorb the new infinities. In SUl ® SUr one only needs counter-terms of a chiral order 
higher than the terms involved in the loops. In Ul®Ur we find a combination of both previous cases. 
It is non-renormalizable and there is a O(p ) term (|i~9|), included to reproduce the ?7a(1) anomaly, 
which at one-loop induces a mixture of chiral orders in the divergent parts, as can be seen from the 
heat-kernel expressions: the divergences are proportional to a 2 (SO) and a = <7( ) +f (2) decomposes 
in two pieces, O(p ) and 0(p 2 ), respectively. (There are divergences proportional to R 2 too 
), but the curvature R associated to the connexion (p§]) does not get any O(p ) contribution). 
There is a lot of freedom in deciding the prescriptions of what removes which divergences, all 
of them equally acceptable from the point of view of rendering the final result finite. They are 
not completely arbitrary, though, since the nesting of divergences when higher loops are considered 
imposes some constraints among the results at different orders, of the Gell-Mann and Low type in 
QED |[10ft. Furthermore, some of them appear more natural than others. 
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Let us analyse the question of the renormalization of the pseudo-goldstone boson masses, in- 
duced by a quark-mass term. Let us first disregard the £7a(1) anomaly by freezing the functions 
of X to their constant values at X = and, for simplicity, consider the symmetric case where the 
three quark species are degenerate in mass, and switch the external sources off: x is a constant 
that multiplies the unit matrix; at tree-level x gives the nonet mass. Now, let us write all the terms 
quadratic in the fields with at most two derivatives, having added the one- loop divergences to the 
tree- level result (prior to renormalization). There are contributions from the operators Ex, Ei in 



(33) for the tree-level parts, whereas the divergent parts come with the operators O4, O5, Oq, Os 



Ox7 and Oig, and can be read off from fl35|) and fl43|). It yields, 

1 - Q^vf • - i X 7? • tt) + ^r,°dV - ~ (l - 2n,-^nA 0O ) x(v°) 2 , (46) 

where Aqo is the ultraviolet divergent amount, that in dimensional regularization is essentially -jyr^, 
(see 36). The difference between the singlet and the octet is apparent. The piece d^rpd^rp does not 
get any divergent contribution while the octet counterpart d^n -d^n does, and exactly by the same 
amount as the mass term tt ■ 7? does as well. In the octet sector, one can pull out the common factor 
from the entire kinetic term and render a finite result by a field redefinition tt — ► (1 — hnijtX^ir. 
There is no infinity left over that could require mass renormalization. 

In chiral perturbation theory it is often simpler to talk about a renormalization of an 0(p 4 ) 
operator rather than a wave-function renormalization, but in our case this is what it corresponds 
to, and this precision is required to qualify the mass renormalization issue. The same result also 
holds in SU L g> SU R . 

For the rj°, though, the divergent contributions to d^rpd^rp - which are none -, and to (rp) 2 
are different, and therefore the mass gets necessarily renormalized by an infinite amount. This is a 
remarkable difference between SUl <8> SUr and Ul <8> Ur. 

In both cases, of course, the divergences disappear if the quark-mass is turned off x ~ ¥ 0> 
which reflects the fact that the spontaneously broken symmetry is built-in, loop by loop, in the 
quantum theory and prevents the goldstone particles from acquiring a mass. When the quark- 
mass is turned on, it is not true that the symmetry structure prevents the rp mass from being 
infinitely renormalized, as happens with the octet mass. The difference can be rooted to the terms 
in the lagrangian that are responsible for the wave-function renormalization. To one-loop, this 
only includes the terms that are quartic in the fields with two derivatives, which are obtained by 
expanding the operator (C^C^); by contracting the two fields that carry no derivatives a tadpole 
diagram is generated and its divergence multiplies d^rr ■ O^tt. The interesting point is that in the 
chiral lagrangian such a term involving the rp field, like (rp^d^rpd^rp or 7? • itd^rpd^rp , does not 
exist at all, as can be immediately seen by making the invariant decomposition of the matrix U in 

U = e V 1 U s , where U s contains only the octet fields and has det U s = 1. then reads 



ni f 



and (C^C M ) 

2 



(CCV) = -—dPrPdrf + (U S %U S U$d»U a ), (47) 

and no crossed singlet-octet term survive for (U^d^U s ) = 0. (C^C^) provides the rp with a 
kinetic term and nothing else, and no term can generate a one-loop wave-function divergence for 
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it. Whereas for the octet, one learns from 

UldfJJs = i jd^-rr + ^[vr^vr] - i --^ [it, [ir,d^}] + ... 

that 

- ^(U$d»U s U%U S ) = l -d^ + ^(M^IMN) + ... . (48) 

It is the second term in (||^) which is responsible for the one-loop wave-function renormalization. 
The structure of commutators makes one realize again that such terms vanish for the singlet. 

This same piece generates also a divergence to the mass term. It gets additional 0(p 2 ) con- 
tributions from operator E2 in (|33|) which come from tadpole diagrams too that arise from four 
meson interaction terms ~ -p(<I> 4 ). It is a well-know result that a lagrangian for scalar fields with 
no derivative couplings other than the kinetic energy with a quartic interaction has an effective 
action that needs a wave-function renormalization that starts at two loops At one-loop it 

requires mass renormalization and this is what we find for the rp field. (Vertices with more than 
four fields from E2 at one- loop do not participate in the renormalization of the kinetic terms). 

We see that the singlet vs. octet difference in mass renormalization is imbued by the structure 
of the symmetry group. 

If one includes the Ua (l)-anomaly effects, both the octet and the singlet get their masses 
renormalized by an infinite amount. 

In this section we have presented the complete one-loop calculation of the divergent part of 
the effective action. It is a calculation to all orders in 1/N C , in the sense that the functions of X 
have been kept generic through the end. There are many unknown parameters in this approach 
(potentially, all the coefficients of the functions of X) and without any further restrictions we would 
not know how to eliminate any of them. We invoke the 1/N C expansion of QCD and the restrictions 
it imposes on the chiral lagrangian to classify the coefficients according to the maximum 1/N C power 
allowed for each, so as to estimate their size, and with this criterion try to select the fewer terms 
that allegedly bring the main contributions. This will be done in the next section. 

4 The 1/N C expansion. 

The systematic expansion of QCD in powers of 1/N C provides a way of effectively reducing the 
number of constants that intervene in a certain process, once it is decided where to truncate the 
series in 1/N C . If N c is large enough, a few terms will suffice to give a good account of the exact 
result. How much large is large enough is a question hard to assess, for a good reason, that despite 
the simplification the large N c limit entails technically, it remains too difficult to sum the subclass 
of diagrams that survive in the limit and it worsens, if anything, for the sub-leading contributions. 
It is argued that, conceivably, big numerical factors might be accompanying the powers of N c in the 
denominator; in that case a few terms would give accurate predictions even for N c = 3. That would 
explain the remarkable resemblance of many qualitative features and patterns of the leading terms 
to those observed in hadron physics, with N c = 3 ||. At any rate, lacking of any analytical result, 
it is the accuracy of the predictions in explaining the data what could give an ultimate justification 
for the expansion. It is this perspective what launched this project, to set out the basis for the 
systematic study of the predictions that come out from such scenario for soft ir, K, rj, rf so as to 
discern in what processes and to what extent an agreement with experiment holds. 
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On general grounds, Witten || showed that in the N c — > oo limit, if QCD confines it has a 
mesonic spectrum that consists of an infinite number of noninteracting, stable states, with masses 
that have smooth and finite limits. Furthermore, the strong interaction scattering amplitudes are 
given, to lowest order in 1/N C , by sums of tree diagrams with mesons exchanged, which can be 
derived from an effective lagrangian with local vertices and local meson fields. The decay constants 
f's are of order y/W c and a coupling constant for a vertex with k mesons attached to it is of order 

(fc-2) 

N c 2 j i.e., it decreases with the multiplicity of mesons in the vertex, each new meson bringing in 
one additional power of 1/y/ N c . The N c counting rules imply that while large- iV c QCD is a strong 
interacting theory in terms of quarks and gluons, it is equivalent to a weakly interacting theory of 
mesons. The higher order corrections in the 1/N C expansion, in addition to new couplings in the 
effective theory, also include the loop diagrams of mesons, which as in any quantum field theory 
reestablish the unitarity constraints on the amplitudes (cuts, discontinuities across, etc.). Each 
loop of mesons, in the effective theory, contributes an extra power of 1/N C . 

In addition to the mesons there are infinitely many glueball states, which at N c — ► oo are stable 
and noninteracting. The amplitude for a glueball to mix with a meson is of order 0(l/s/Nc), and 
the vertices in which they are involved are even more suppressed than the meson vertices: one 
power of 1/N C for each glueball. 

However, only ir, K, r], rf remain massless in the chiral limit and N c — > oo, their masses being 
precluded by Goldstone's theorem. The rest of excited mesons and glueballs remain massive, with 
a typical hadronic mass of about 1 GeV. They decouple from the soft processes that involve the 
goldstone bosons by virtue of their masses and they are integrated out from the effective theory. 
The baryons have masses much higher for large N c , for they are known to grow like N c ||. 

Within the large- iV c , the Ua(X) anomaly effects can be accommodated in a rather natural way 
in the framework of the chiral lagrangian. Actually, the identification of each ingredient that has 
been taken into account in writing down the effective theory is clearcut: the constraints imposed on 
the interactions among goldstone bosons are contained in the operators that involve the U matrices 
(but not logC/); the quark masses enter through the terms in \\ and the UaO-) enters through 
the functions of X. One can talk of switching off the J7^(l) anomaly by freezing the functions 
to constants, in a similar way as one can take the chiral limit by sending the quark masses to 
zero. In this language, one can say that the rf has two ways of manifesting itself: either as a 
goldstone boson or as the argument of the functions of X, breaking chiral symmetry as dictated 
by the ?7a(1) anomaly. When it manifests as a goldstone boson, its couplings follow the rules of 
the meson couplings. This is its the mesonic part, associated to the content in quark degrees of 
freedom. Its other presence in the chiral lagrangian, imposed by [/^(l), involves couplings that are 

3 

more suppressed, 1/N£ per rf in a vertex. This is associated to the special gluonic content of the 
rf , put forward by the anomaly. Although there would be no such a thing as the rf in a world 
without quarks - nor chiral symmetry -, in the chiral limit the rf mass is S 



m 2 v o 



Ua(1) 




no quarks 

+ 0(— o' 



where E is the vacuum energy in a world without quarks and with a coupling to the topological 
charge Tr G^G^ in the lagrangian, as in (||). The fixed proportion of rf and 6 that appear in the 
combination X = i ^p l rf + 6 actually relates glueball and rf anomalous couplings to operators 

3 

involving goldstone bosons. The vertex suppression of 1/N£ per rf is a combination of 1/N C - 
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glueball and l/v^A^-meson suppression, the former originated in the anomaly equation, the second 
carried by the / factor that divides the rf field. 

In order to obtain the 1/N C power counting of the sub-leading pieces too, one might proceed 
by comparison of Green's functions, as evaluated from the chiral lagrangian and from the diagrams 
in QCD. In the effective theory, given an operator multiplied by a function G(X), the 1/N C power 
counting can be established on the basis of two distinguishing features: the number of traces over 
flavour indices (#(t/)), and the number of powers of the source 6(x), (#(#)). As a rule of thumb, 
the leading dependence on 1/N C of the various couplings follows from the simple prescription [jOI]: 

G(X) = g(X/N c ), (49) 

where g is a function whose expansion in powers of X/N c has coefficients of order 1. The origin 
of each factor can be easily traced: in relation to the vacuum energy which is 0(N C ) 2 , there is 
one power suppression of 1/N C for each flavour trace and one for each source - recall that in QCD 
the sources for 1/N C suppression are the loops of quarks and the non-planarity of the diagrams. 
Each trace taken over flavour indices amounts to a sum over quark flavours, which in turn can 
arise only in a quark loop in QCD and adds a factor of 1/N C . As for the source 9(x), it couples 
to the topological charge in (|6|) with strength 1/N C ; each derivative of the generating functional 
with respect to 9{x) will bring one power of 1/N C to the Green's function. In the effective theory 
this is achieved by pulling out an explicit power of 1/N C for each 9{x). Finally, it was already 
mentioned that the ubiquitous factor of / count as y/N c . In particular, the rf that appears in X 
is always suppressed by a factor 1/f (|i~7|). There are no additional powers of 1/N C for the leading 
contributions: once all the previous factors of 1/N C have been pulled out, also from X as X/N c , 
the expansion of g in powers of X/N c has coefficients that are order 1. 

Recall at this point that the 1/N C counting should be done in a chiral lagrangian with a generic 
number n\ of light flavours. This is because of the n;-dependent factorization relations already 
mentioned in section 3, that give linear relations among the a priori independent operators for 
particular values of n^. The mismatch of powers of 1/N C and the departure from the general rule 
(|49|) are avoided by allowing for a generic n/. The correct counting is thus obtained for the functions 
Pi(X) in p3),(B3). The implications for the Li(X) can be read off from (E 



Furthermore, as pointed out in [24], it is the Ul<S>Uh lagrangian that provides the suitable basis 
to establish the correct 1/N C power counting of the constants that involve the nonet of mesons. 
The rj' in the large-iV c limit gives a contribution to the large-distance behaviour of the Green's 
functions that should not be overlooked; in the limit its properties are the same as for the rest of 
goldstone bosons in the nonet. In next to leading order, a topological mass term appears for the 
rj', but it is 0(1/N C ) and is treated perturbatively. In counting powers of 1/N C , the rj' cannot be 
integrated out from the nonet, for it is when N c — > oo, when the 1/N C expansion is more sensible, 
that the rj becomes massless and does not decouple. 

Expanding the Wj's functions in ( |i~9| ) in power series in X, 

W k (X) = W k0 + W k2 X 2 + W M X A + ... = £ (v kQ + v k2 X 2 + v M X 4 + ...) , (50) 
for k = 0, 1, 2, 4, 5, 6, and for W 3 

W 3 (X) = W 31 X + W 33 X 3 + ... = -*£ (v 31 X + v 33 X 3 + ...) , (51) 
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and following rule 



we find 

W 00 = 0(N 2 C 

Wxo = W 2 o 

W 31 = O(l), 



4 

^50 = O(l) 



T-r, 



60 



0(1) 



For the coefficients of the 0(p 4 ) operators that do not involve derivatives of 6 [24|, for 
we find, 



£i(0) ,L 2 (0) ,L 3 (0) ,L 5 (0) ,L 8 (0) ,L 9 (0) ,L 10 (0) ,L n (0) ,L 12 (0) , 

Li 3 (0) ,L 14 (0) ,Li 5 (0) ,L 16 (0) = 0(JV c ) 
2Li(0)-L 2 (0) ,Li 3 (0)+8Li(0) , L 14 (0) - 2^(0) ,L 4 (0) ,L 6 (0) , 

L 7 (0) ,L 18 (0) ,L 19 (0) ,L 20 (0) = O(l) 

L 15 (0) - 2L 14 (0) , 2L 16 (0) - L 14 (0) , L 17 (0) = 0(l/iV c ). 



(52) 



For the parity-odd terms, the first contribution is given by the linear term of the X expansion, 

L 21 (0),L' 25 (0),L 28 (0) = 0(1) 
L 22 (0),L 23 (0),L 26 (0),L 27 (0) = 0(1/N C ) 

L' 24 (0) = 0(1/N 2 ). 

(53) 



2 

m n 



As mentioned, the r]° gets a contribution to its mass that is 0(1/N C ), in the notation of ([5l] ) 
= —nivo2- Counting m 2 as two chiral powers 0(p 2 ), in the multiple expansion we shall 

count 1/Nc also as 0(p 2 ) [[0|. However, to be fully consistent with it, if terms 0(p 4 ) x j^- are 
kept, then the chiral 0(p G ) order should be also included. This would require a two-loop chiral 
perturbation theory calculation which is far beyond the scope of this article. 

Following these criteria and using (|50|), ( |5l| ) and (|52|), we expand the Bi and Oj functions from 
Appendix D, first in powers of X and then in powers of keeping corrections up to jj- for the 
0(p 4 ) terms, up to for the 0(p 2 ) terms and up to -tj-j for the O(p ) one. Notice that none 

iVc -**c 

of the rf fields that appeared through the X has survived, which means that all the terms in the 
lagrangian starting with more that four fields are eliminated. This kind of terms would be only 
required for calculating processes that are very difficult to measure experimentally. 
To this order, only two fVs survive: 



02 



\ V02 + n\ v 2 V31 + 0(j± 



(54) 



The Bi's for i=0 to 20 are the same as in (f43|) except for two new contributions: 



5s 



N c 2 



B 



12 



(55) 
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The rest of coefficients either vanish exactly or do not contribute to this order in the expansion. 
There are also some contributions that are proportional to 9. 

The use of the equations of motion does not ruin the 1 /N c power counting in (|49| ) , although it 
introduces additional products of traces and powers of derivatives of 9 in (^). The very structure of 
( |23|) complies with (|49| ) since the functions of X that multiply the various operators carry their own 
1/N C suppression required by ©: for the identity, jjji = g(X/N c ); (C"C M }, (C") C M and D»d 
are multiplied by ^ = ^- g(X/N c ); N and M do not get any suppression since both ^ and 

are g(X/N c ); however, (M) and (iV) appear multiplied by and |^ which are g(X/N c ); 

= g(X/N c ) multiplies D^D^ whereas ^^^e = g(X/N c ). Here 5 denotes, in each case, 
some function that has iV c -independent Taylor expansion coefficients. 

Finally, let us comment on the regeneration through quantum corrections of a term 
(U^ D ^U) {U^ D^U) , which had been removed from the tree level effective lagrangian in (|l9|). It 
is a confirmation that the effects of meson loops bring to the effective action contributions that are 
l/N c suppressed in relation to the leading tree level. It is only at N c = 00 that ir, K, 77, rf form a 
nonet. When sub- leading corrections in 1/N C are taken into account, the enlarged symmetry with 
respect to Ul <g> Ur no longer holds. The first instance is the 0(1/N C ) mass piece exclusive for the 
r/o- The reappearance of that term is another example: the rp (singlet) and the ir (octet) fields are 
normalized differently by sub-leading contributions. 



5 Conclusion and outlook. 



In this article the effective theory developed in [If] for the strong interactions at low energies 
among the lightest pseudoscalars is extended to include the rf particle. The approach that has 
been adopted exploits the fact that, as the number of colours N c grows bigger, the mass difference 
(topological mass) between the rf and the octet vanishes as 1/N C , and the Ua(X) puzzle can be 
treated as a series in inverse powers on N c . We exploit the possibility that a good description 
could emerge by taking the nonet of soft pseudoscalars as the goldstone bosons of the spontaneous 
breaking of Ul £3 Ur — > Uy, thus benefitting from the perks that such theories feature, in terms 
of constraints on the form of the couplings and relations among the coupling constants. The 
departures from the results in the real world are dealt with as corrections in powers of the quark 
masses and 1/N C . 

The most general effective action is given that includes terms up to 0(p 4 ), with quantum 
corrections included to one-loop. The Ua(1) anomaly is conveniently incoporated. 

This is the first step of a project towards a systematic study to spell out the predictions that 
could emerge from such scenario. It will shed light on what processes can be accomodated within 
the approach. 
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7 Appendix A: The Heat-Kernel Technique. 



This appendix is a reminder of some results that have been used in the text, which are obtained 
with the help the heat-kernel technique. It provides a convenient way to evaluate the one-loop 
effective action. 

Let H be an operator acting on a Hilbert space. The heat-kernel two-point function associated 
to H is defined function of the parameter r, 

K(x,y,T) = (x\e- iT& \y)e(T). 

It verifies the equation 

i^-K(x,y;r) = J d D z (x\H\z)K(z,y;r) + i5(r)5 D (x - y), 
with the boundary condition, 

lim K(x, y; r) = 5 D (x - y), 

which follows from its definition. We do not specify the dimensionality of space-time and allow for 
a generic D. Often, H is a local operator, i.e., (x\0\y) = O x 5 D (x — y), where O x is a differential 
operator with a finite number of terms. We shall limit ourselves to this case. Then, K(x, y; r) is a 
Green's function that verifies a partial differential equation of the Schrodinger type, 

d 

i 0^ K ( x ^ y> r ) = °x K (x, y; t) + i5(t)5 d (x - y). 

When the operator H is elliptic the definition of K differs somewhat from the one given above. 
In that case there is no need for a factor of i in the exponent, since all the eigenvalues of H are non- 
negative. There, the equation is a heat-transport-like equation, from which the technique shares 
its name. In the present case we are dealing with operators of the sort d^d^, which in Minkowski 
space is hyperbolic. It is convenient to pick out the mass term from H, if the theory is massive, or, 
else, to introduce an infrared regulator by adding a constant term M 2 to H, H + M 2 . For small 
values of r, K(x, y; r) admits an asymptotic expansion of the form, 

1 r - ) 2 

K(x,y;r)~i ^ e "^+^ £ h n ( x , y )(i T )« . ( 56 ) 

(4vrzr) 2 n=0 



The functions h n (x,y) are known as the Seeley-DeWitt coefficients [28] 
Given that 



1 (x- y ) 2 n 
K(x, y; t —> 0) ~ % 5- e 4iT — ► 5 (x — y), 



(A'Kir) 

the boundary condition translates into 

ho(x, x) = 1. 

The computation of the one-loop effective action using the background field method entails the 
evaluation of 

logdetCff) = Tr log If. 
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Tr stands for the trace over space-time as well as all the internal indices. This can be written in 
terms of the kernel K(x, y; r) if one uses the following integral representation for the logarithm 



logH- logtfo = -|°°y (e~ iTii - e~ iTAo ) . 



Apart from an unessential H- independent, divergent constant C 

/poo J T 
d D x J — tr (x\e- lTH \x) + C. (57) 

where tr stands for the trace over internal indices only. 

This expansion will be used to depict the divergent short-distance contributions so as to be able 
to subtract them away. Given r, only those points separated an interval (x — y) 2 that is less or of the 
order of r give non-negligible (non-oscillatory) values for K(x, y; r). The short distance contribution 
to the effective action is thus contained in the integral about r = 0; it will not be affected by M, 
which only modifies the contribution from large values of r (large wavelength modes). 

The particular form of the operator H which corresponds to our computation is 

H = + a , 

where d M = d^+uj^x). A naive application of eq. (|57| ) could lead to infrared divergences associated 
to the large r integration region. For that reason we have changed it to H + M 2 = d^cP + a + M 2 . 
Each term in the Seeley-DeWitt expansion yields a finite contribution, which give 



1 ^ f D \ 

i tr(x\ log (H + M 2 ) \x) = E( M ') f [ n ~ y J (hnfa x ))> ( 58 ) 

except for an unessential additive constant. The Seeley-DeWitt coefficients are known in this 

case, 

h (x,x) = I, 
hi(x,x) = —a, 

h 2 (x,x) = loZ + ^R^Rr + ^tPa. (59) 

where R^ v = d^ujy - d^u^ + [ k>^,k>i/]< Higher order h n (x,x) may be found in |28| . 

The singularities that arise in the calculation have been regulated by analytic continuation over 
the D-complex plane. In a four-dimensional theory the ultraviolet divergences appear as poles 
about D = 4, and get contributions from ho, h\ and hi, as T T (l — ^) an d r (2 — y 

respectively. Retaining in (pq) these three terms only, it reads 



i tr( 



<x|log(H + ^)|x>| = _(_ + log4 J r- 7+ l <-(.+ M 2 f + -R^ 



(4vr) 2 V^>-4 J x 2 V 7 12 

+ Wf { (l M4 " H " T2 R ^ U > + ° {D ~ 4) ' (60) 



^In general cr{x) and u)(x) can be matrices. In that case the coefficients h n are matrices as well and a trace over 
these internal indices is also understood in (158). 
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Notice that in fl60| ) it is the combination a + M the one that multiplies the pole 7734, i-e., 
it is the whole operator that adds to dpd^ in H + M 2 . If a has its own mass terms - the light 
pseudoscalar masses in our case - there is no need to introduce any new infrared regulator M. We 
see that only is involved in the residue of the pole ■ 

Therefore, a finite expression is obtained by subtracting 



+ log47r-7 + l ( h 2 (x,x) ) 



(4vr) 2 V D - 4 

to the effective action. This copes with the ultraviolet divergences and it is the procedure used in 
the text, fi is a parameter with mass units. The total derivative (dpd^a) has been discarded. 

8 Appendix B: U(3). 

This appendix is devoted to present some of the properties of the U(3) group which have been used 
in the text. Whenever the generalization is straightforward we give the result for U(n). 



The explicit form of the U(3) hermitian generators , A M = At , \i = 0, 1, 2, 



is 








i | , (61) 




/ 1 












—i 











i 



















(° 



















, A 7 = 
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1 



~T 1 
^ 3 » -2 



For U(n) there are n 2 such matrices, and \i = 0, 1, ...,n 2 — 1. Also Ao = 
identity. These matrices obey, in general, the following basic trace properties: 



>?m 5, 



no 



(A M A V ) 
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~I, with I the n x n 



(62) 



Tr (A M ) = (A^) 
The product of two matrices verifies: 

[A^, X v ] = 2if /Jl i,p\p , {A/i, A^} = Id^pXp , 

A^jAjy = {dpyp ~\- if pup) Ap . 

For U(3), the non-zero entries of the antisymmetric fp Up and symmetric dp Vp constants are 



(63) 
(64) 



/l23 = 1 , /458 
/l47 = — /l56 = /246 = /257 = /345 



7678 — — , 



/367 



1 

2 ' 



(65) 
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dll8 = d,228 = d%38 = —dg88 = TJTJ , 
^146 = ^157 = —^247 = ^256 = ^344 = ^355 = "^366 = _ ^377 = \ i 
(i448 = C?558 = C?668 = ^778 = — 2 "v 7 3 ' 

These tensors are not independent as they satisfy the following relations (repeated indices are 
summed up): 

d^v = n^/2n 5p , 
d^ v \dpv\ = n (6fj, p + S^oSpo) , 

fpurfxpr fpXrfpur ~H fpprfuXr = > 
ffiurd\p T + ffiXrdpur ~\~ j \iprd u \ T = , 



ffivcrfpra — dpp a d T y a dpT-Q-d^pu . (67) 



The identity 



(A a ) afe (X a ) cd = 2<5 a Ac , (68) 
has been extensively used. It yields the properties 

X a X a = 2nl, 
(X a AX a B) = 2(A)(B) , 
(X a A)(X a B) = 2{AB) . (69) 

9 Appendix C: Operators with derivatives of 9 and some 
C-violating operators. 

In this appendix to give list of the operators that involve derivatives of 6, that we have not included 
in the text (see (|32|)). 

031 = Bpfi {C»C V C V ), 

32 = D^O (CWC), 

033 = Dfft (C»C)(C U ), 

034 = Dp§ (cn(cn(c„), 

035 = D^O D»§ {CC V ), 

36 = Dp§ dJ (c^cn, 

037 = D„0 D»0 (C V )(C V ), 

38 = Dp9 Dj (C»)(C), 
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o 39 
O40 
041 

A2 

043 
044 
045 

O m 
O47 
48 
O49 
O50 
O51 
O52 
O53 
O54 

O55 

057 



£> M D*d D V 6(C U ), 

D^e D^e dJ d u §, 

iD^D^O (C U C U ), 

iic^D^e d v dJ, 
iD^B D^O d u dJ, 

D^D^O D U D U §, 
D^e {CM), 

d^o (C)(m), 
iD^e {cm), 

iD tl 9{C){N), 

d^o d^o (M), 
iD^e (n), 

iD^D^O (M), 
D^D^O (N), 

D,6 (C v )(Ff 



The left-over independent operators are C-violating and are the following: 



(C lt )((F£' + tfFg'u)c v 
{F^F Lllu -F^F R(lv ), 



(70) 



(71) 



-pupa 



{F£ v + Fg){CP) D a 



10 Appendix D: The one-loop renormalization functions 



In this appendix, we collect the functions fVs an d B^s, as defined in equation (35) , which have 
been computed with the heat- kernel method. The results have been obtained with the help of the 
program Mathematica. 

The 0(p ) lagrangian is renormalized with Qq, whereas the rest of the Qj's renormalize the 
0(p 2 ). For technical reasons, it was convenient to keep one more operator at 0(p 2 ) than those 
given in (|33"|), E 7 = iD^D^O. By integration by parts, this can be written in terms of the other Ei, 



fl 7 E 7 = -in' 7 E 5 - iQ! 7 E & 
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which should be taken into account and added to and Qq. 
For the sake of brevity, we use the following short notation: 

U -Wi ,_w< _ w>> 



and, similarly, to include as many primes as necessary. Let us emphasize that ui^ is not the derivative 

f 
4 



r'2 

of Ui, and so forth. Also, recall that W^i(O) = ^(O) = v so that u>i(0) = oj 2 (0) = 1. 
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iZ o - -g- ^0 — g- u o u o 1 "32" 1 

O — n ^ , , n ^ , ," , ," _l_ n ^ , ," , / 2 , n f , / , / , ,// n t , J , ,/ 3 

III — W W l 7f ^0 W l + — W W l + — ^0 1 1 ~ ~[Q 1 ' 

_ 1 // 77-7 »» »» . 77-/ a f \ 1 1 Tli 1 1 

\l 2 = --U> U 2 — LO i0 2 + ^0 ^3 + 4 ^0 w l ^2 g- ^0 ^1 ^2 

2 2 
^/ I I n / 11 1 1 • 111 • ^/ // / 

+ -g" w w l ^2 + -g- ^0 w l w 2 - ^0 w l ^3 - *"g- w ^1 ^3 

4 3 
77*7 » ill • ^7 / / 2 

- ^0 w l w 2 "^Tg w l w 3, 
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+ 



o ■ n l 11 1 1 // n l 11 11 . ■ n l 11 1 . ■ n l 1 1 1 

"3 = -«T- W w 2 - T w w 3 J" ^0 w 3 +*-^^0 w l ^2+^-j-^o w l ^2 



I - - — I W 0^ W 3 + — O; CJ 3 + — W - — CJ w 3 



• I 12 

- uj uji io 2 , 
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11 Appendix E: The Discrete Symmetries: C, P and T. 

For the sake of completeness we give in this appendix the transformation laws for the fields and the 
sources that appear in the article under the discrete symmetries, C, P and T. In order to specify 
them one needs to first define how they act on the space-time coordinates x^ = (t, x). 

The Parity (P) operation transforms x — > — x while leaving the time component unchanged. 

Using the Minkowski space notation, x M — ► x^ = p^ u x u , where pP v = diagonal (1, —1, —1, —1) is a 
matrix. 

Time-Reversal (T) reverses the flow of the time-component t —t while leaving the space 

components unchanged, x^ —> = t^x u , where t% = diagonal (—1, 1, 1, 1). 

Charge-Conjugation (C) does not act on space time-indices, it interchanges the role of particles 
and anti-particles. 

In Quantum Mechanics they arc implemented with operators acting on a Hilbert space that are 
unitary for C, P; and anti-unitary for T. 

Acting on the (Dirac) quark fields q a (x), where a labels any colour or flavour index, they read 

qa (x)^ q P( x ) = ZcCqa T (x), C 7 JC- 1 = - 7/1 ; 
-^q ( a P) (x) = Z P Pq a (x p ), Vi\V- x = lv ;, 

-Z+q<P(x) = Z T Tq a (x t ), TtfT-^^, (74) 

The £cs £p> £t are arbitrary phase factors, \£,c\ 2 = |£p| 2 = |£t| 2 = 1- The matrices C, V, T act on 
Dirac indices only. In the Dirac representation C = «7°7 2 , V = 7 . Once C, V are fixed, T = —i^C. 
They verify C" 1 = C 1 " = C T = -C, and T -1 = = -T T = T. Acting on 75 they yield 

C7IC- 1 = 75, 7°757° = "75, r^T = 75. 
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For q a (x), 

C 

CpQa(Xp)j 



q a (x) — > -CclI( x )C \ 
p 



&g (zt)T, (75) 

The phase factors £q, £p, £y shall be omitted henceforth. The quark bilinears q a (x)Tqb(x) transform 
as 

q b (x)[T] c q a (x), [T] c = (c^Tcf , 

q a (x p )[T] P q b (x p ), [T] P = 7 r 7 , 

g a (x t )[T]T56(xt). [r] T =Tr*T. (76) 

The star (*) in the last line of ([7(]) denotes complex conjugation. There is a minus sign in the 
bilinear transformed under C which comes from the anti-commutation of two quark fields. 

[I\c = I, [I] P = I, [I\t = I, 

[ilb]c = «75, [ilb]p = -^75, [*7s]t = -i 7 5, 

[f] c = -y*. [f] p = p ^'i [rf] T = -t^W% 

[7 M 7s]c = [7 M 7s], [7 m 7s]p = -J^b^Ts], [7 M 7s]r = -^[V'ts]- 
For the gluon field (hermitian) matrix G^{x) in colour-space, 

^-G^ T (x), ^p;,G"'( Ip ), £-^G*»'(s t )- 

It is easy to verify that the QCD action is invariant under C, P and T. 

For the topological charge density Q{x) Tr c G^(x)G pa (x), which is also real, 

%Q{x), £ det(p^,) Q{x p ) = -Q(x p ), ^ det(^,) Q(x t ) = -Q{x t ). 

So far for operators involving the dynamical fields. 

The (hermitian) operator iq a ( x )l5qb( x ) has the same quantum as the light pseudoscalar matrix 
& a b(x), and, since the vacuum is invariant under C, P, and T the transformation laws of the latter 
are taken from those of the former. Let us write them down for the simpler case of Ul{2) <g) Ur(2), 
where 

' TT 



* = f- _zm±zz I • ( 77 ) 



Under the discrete symmetries 



n V2 



$ T , (78) 



and, similarly, 



*(x p ), -^-$(z t ), (79) 
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which generalize immediately for Ul(jh) ® Uji(ni). It translates into 



U(x) 




= f/ T (x) 


U{x) 


^ C/( p )(x) 




U{x) 


C/ (T) (x) 




in f2|). 







(80) 



As for the external sources, we shall chose their transformation so as to leave the action invariant. 
The real source 9(x), that couples to the topogical charge Q(x), transforms as: 

£ e(x), £ -e(xp), ^ -e(xt). 

In the text, the combination 9 = iO appeared in a natural way. It transforms accordingly, with an 
extra minus sign for the T transformation because it anti-commutes with the imaginary number i 
due to its anti- unitary character. The combination X, defined in (|i~7|), transforms as 9 does. 
The (hermitian) source matrices s ,p , a«, in flavour space transform as 

s(x) -> s T (x), s(x p ), s(x t ), 

p(x) -> p T (x), ^ -p(xp), -p(x t ), 

£ -^ T (x), ^p^t/fr,,), - -^^'(x t ), 

The combination % = 2B(s + ip) transforms as the U fields. 

The left and right combinations of the vector and axial sources, = v M — a M and r M = + a^, 
transform as 

£ -r^ T (x), ^^^'(xp), - -^^'(x t ), 

r"(aO £ -/^ T (x), ^p^'(xp), £ -t>"W 

Both the C and the P transformations interchange left and right. 
For the field strengths F£ v , F£ v associated to 1^, r^, 

F^(x)Z-F^ T (x), 4 P y,Ff'( Ip ), ^-fyt^f'te), 

F^(x)Z-F^ T (x), Z P ^,F^\x p ), ^ -t^ vl F^\x t ), 
Finally, for the combination C M = U^D^U, that is anti-hermitian t = — C 1 , 
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